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Abstract 

Coarse-grained phase distributions are introduced that approximate to the 
Susskind-Glogower cosine and sine phase distributions. The integral re- 
lations between the phase distributions and the phase-parametrized field- 
strength distributions observable in balanced homodyning are derived and 
the integral kernels are analyzed. It is shown that the phase distributions 
can be directly sampled from the field-strength distributions which offers the 
possibility of measuring the Susskind-Glogower cosine and sine phase distri- 
butions with sufficiently well accuracy. Numerical simulations are performed 
to demonstrate the applicability of the method. 



1 Introduction 



Although the problem of defining and measuring quantum mechanical phases 
of radiation-field modes has been discussed for a long time, there has been 
no unified approach to the phase problem so far [l| • The reason is that there 
is no "good" (i.e., self-adjoint) phase operator in the Hilbert space. The 
difficulties obviously arise from the boundedness of the eigenvalue spectrum 
of the photon-number operator, which is the canonically conjugate of the 
phase operator. 

One way to overcome the difficulties was proposed by Susskind and Gl- 
ogower |2j who used the non-Hermitian exponential phase operator and its 
adjoint in order to define two Hermitian operators. From the analogy be- 
tween them and classical trigonometric functions the two operators are also 
referred to as cosine and sine phase operators. In particular in the classical 
limit they exactly correspond to the cosine and the sine of the phase. Since 
the cosine and sine phase operators are self-adjoint, their eigenstates can be 
used to define proper probability distributions for observing the cosine and 
sine phases. However, these distributions cannot be measured, in general, 
simultaneously and do therefore not uniquely characterize the phase of the 
quantum state of a mode. This is obviously the price to pay for introduc- 
ing self-adjoint operators and keeping the concept of quantum mechanical 
probabilities. 

In the classical limit the difference between the cosine and sine phase 
distributions vanishes. In particular, in classical optics the phase difference 
between two radiation-field modes can always be determined by simultane- 
ously measuring the cosine and sine of the phase difference in interference 
experiments. Inspired by such kinds of measurements, Noh, Fougeres, and 
Mandel || introduced an operational approach to the cosine and sine phase 
operators. As expected, for classical fields the corresponding phase distri- 
butions agree with the quantum mechanical Susskind-Glogower (SG) phase 
distributions, but for quantum fields they may be quite different from each 
other. 

A powerful interferometric method for determining phase-sensitive prop- 
erties of optical fields has been balanced homodyne detection. It is well 
known that from the data recorded in a succession of measurements the 
quantum state of a signal mode can be obtained, as was first demonstrated 
experimentally by Smithey, Beck, Faridani, and Raymer using optical ho- 
modyne tomography (|J. Knowing the quantum state (in terms of, e.g., 
the density matrix in a chosen basis), the quantum phase statistics can be 
calculated ||. In the present paper we show that the SG cosine and sine 
phase distributions can be measured more directly by sampling them from 
the recorded data |§. 

The approach is based on the introduction of parametrized cosine and 
sine phase distributions that are defined on the basis of appropriately coarse- 
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grained SG cosine and sine quantum states. So, the parametrized distribu- 
tions can be regarded as smoothed SG cosine and sine phase distributions 
that reflect the exact statistics in a coarse-graining approximation. The ac- 
curacy is determined by the parameter that defines the phase coarse-graining 
interval. In particular, the parametrized distributions tend to the exact ones 
when the coarse-graining parameter approaches zero. 

The advantage of the method is that for any non-zero value of the param- 
eter the coarse-grained sine and cosine phase distributions can be directly 
sampled from the recorded data in balanced homodyning. Hence, choos- 
ing the coarse-graining interval small enough, the exact distributions can be 
measured with sufficiently well accuracy. Since in balanced homodyning the 
observed difference-count statistics is (for chosen phase parameters of the ap- 
paratus) a scaled field strength distribution of the signal mode, the sampling 
function needed is given by the kernel in the integral relation between the 
parametrized cosine and sine phase distributions and the phase-parametrized 
field-strength distributions. 

The paper is organized as follows. In Sec. ^| coarse-grained cosine and 
sine phase distributions are introduced. Section ^| is devoted to the relations 
of the phase distributions to the field-strength distributions and the calcula- 
tion of the cosine and sine phase sampling functions. In the Sec. |] results of 
computer simulations of measurements of the cosine and sine phase distribu- 
tions for a squeezed vacuum state are presented, and a summary and some 
concluding remarks are given in Sec. ||. 

2 Coarse-grained cosine and sine phase dis- 
tributions 

Using the exponential phase operator 



n=0 

(n \n) =n\n), h = a) a, [a, at] = 1), Susskind and Glogower p| introduced the 
Hermitian cosine and sine operators 



oo 



E =J2\n)(n + l 



(1) 



C 



i(£+£t) 



(2) 



and 



S = -±i(E-fr), 



(3) 



respectively, which satisfy the eigenvalue equations 




(4) 



(0 < < tt) and 




(5) 
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(— 7r/2 < 0< vr/2). The SG cosine and sine phase states | cos0) and | sin0), 
respectively, form orthonormal and complete sets in the Hilbert space, that 
is to say, 

(cos 0| cos 0') = 5(0-0'), (6) 



d0 | cos 0) (cos 0| = I (7) 



and 



(sin 0| sin 0') = 5(0-0'), (8) 

C7T/2 

-vr/2 

In order to give a unified approach to the states, let us consider the ip- 
parametrized Hermitian operator 



/■tt/2 

/ d0 I sin 0) (sin 0| = I. (9) 

J-tt/2 



C{$) = \{Ee~^ + £t e #). (io) 
Recalling Eq. ([[]), C(ip) is easily proved to satisfy the eigenvalue equation 

C7(V0|$,V) = cos$|$,V>, ( n ) 



where 



W- Ve im/, sin[(n + l)$] |n). (12) 
V vr n 

Note that | — $,-0) = -|$,^) and |$ + 7r,-0) = |$,-0 + 71 ")- For chosen -0 the 
states |$, i/j), 0<$<7r, form an orthonormal complete basis. The SG cosine 
and sine phase states in Eqs. (f|) and @, respectively, can be obtained by 
appropriately specifying the states \Q,ip), 

|cos0) = |$ = 0,^ = O), (13) 

|sin0)= $ = ±vr-0,^ = ±7r^ . (14) 
With regard to measurements, we now introduce coarse-grained states as 

|$,V,e> = -^ J d$'|$',^>. (15) 

$-e/2 

Using Eq. fll2|), the states |$,-0,e) can be given by 

|$,<0,e) = W — ^e m ^sin[(n + l)$]sinc[(n + l)e/2] |n) (16) 

' n 

(sine a; = sin x/x). They are normalized to unity, 

($,^,e|$,^,e) = 1, (17) 
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and tend to the exact states \&,?p) as e approaches zero, 

limine) = (18) 

The coarse-grained cosine and sine phase states, respectively, are obtained 
as | cos0, e) = |$ = 0, ^ = 0, e) and | sin</>, e) = |$ = |-7r — 4>, ip = |7r, e) [cf. 
Eqs. flT3|) , (fnp , and (P~5[)l . According to Eq. (|T8|), they approach the exact 
SG cosine and sine phase states in the limit when e— >0. 

The states |$, ■?/>, e) can be used to define parametrized phase distributions 
of a radiation-field mode via their overlaps with the quantum state g of the 
mode, 

p($, e) = iV" 1 ^, e) ($, V, e| g |$, V, e), (19) 
where the normalization factor 

iV(^,e)= f d$($,<M£|$,^,e> (20) 
Jo 

has been introduced. Note that p(Q, if), e) — >p($, </>) = ($, ^I^I^j when e— >■ 
0. In particular, the coarse-grained cosine and sine phase state distributions 
p c (4>, e) and p s (4>, e) are given by 

p c (0,e)=p($ = 0,^ = O,e), (21) 

p s {4>,e) =p($ = l7r-0,^ = |vr,e) (22) 

[Pc(0) = (cos0|^| cos0) = lim e ^o Pc(<f>, e), Ps(</>) = (sin0|f?| sin</>) = lim e ^ 
p s (0, e)]. The smaller the value of e becomes the better the coarse-grained 
distributions approximate to the exact ones. In Figs. 0- |2] plots of p c (4>,e) 
and p s (4>, e) for various values of e are shown for a mode prepared in a coher- 
ent state. We see that with increasing value of the mean number of photons 
the value of e must be decreased in order to obtain the coarse-grained dis- 
tributions comparably close to the exact ones. This is of course a reflection 
of the fact that with decreasing value of e the value of the (effective) cut-off 
photon number in the expansion ( |TE| ) is increased. 

3 Relations to the phase-parametrized field- 
strength distributions 

When we let e = 0, then Eq. ( |TE[ ) reduces, after multiplication by e -1 / 2 , to 
the expansion in the photon- number basis of the exact states \<&,ip). All 
the photon-number states are seen to contribute to the exact states with 
comparable weight, which prevents one from directly sampling the exact 
distributions from the difference-count statistics recorded in balanced homo- 
dyning. The advantage of the coarse-grained distributions is that they can 
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be directly sampled from the recorded data. Since for sufficiently small val- 
ues of the coarse-graining parameter the measured distributions approximate 
to the exact distributions with desired accuracy, the method enables us to 
asymptotically measure the exact distributions. 

As already mentioned, the difference-count statistics recorded in balanced 
homodyning represents the statistics of a scaled field strength of the signal 
mode, 

F(<p) = iFKae-^ + o+e^), (23) 

where the phase (p is determined by the chosen phase parameters of the 
apparatus (see, e.g., ]/]). The desired sampling functions can be therefore 
obtained from the integral relations between the coarse-grained phase distri- 
butions p(<&, ip, e) and the phase-parametrized field-strength distributions of 
the signal field, 

p(J 7 ,<p) = (J 7 ,<p\Q\J 7 ,<p), (24) 
ip) being the eigenvectors of F(ip) (for details, see J7], ||). For this purpose 



we note that the single-mode density operator can be expanded as || [TO], |TT 

rn foo 

Q= dip dJ r p(J r ,(p;s)K(J : ',ip;—s), (25) 

JO J-oo 

where the smeared field-strength distributions p{T , <p;s), s = 1 — 77 -1 , have 
been introduced which are measured in non perfect detection, i.e., when the 
detection efficiency 77 is less than unity (see, e.g., 0). The operator integral 
kernel K{T ', (p; —s) in Eq. ( f2~5"D is given by 



k(F,ip;-s)= / dy\y\ expUy \f&)-F\ -\sy 2 \F\ 2 }. (26) 

Using Eq. (P5|), from Eq. ( |T9"D we easily find that p(<&, if), e) can be related to 

p[T , (p; s) as 

p($, V, e) = N- 1 ^, e) dip dTp{T, <p; s) K e ($, V, J 7 , <p; s), (27) 

J J-oo 

where 

K e ($, V>, T, ^ s) = ($, V, e\ K{T, ^ -s) |$, V, e). (28) 

Equation fl27|) can be regarded as the basic equation for direct sampling 
of the phase distributions ip, e) from the difference-count statistics in 
balanced homodyning, where the integral kernel play the role of the sampling 
function. In order to calculate them, we substitute in Eq. Q2"8D for \&,ip,e) 
the expansion Eq. (|16|), and we derive 

Ke($,ip,F,<p;s) = — $3 53 {fnm(x;s) exp[i(n-m)(<p-ip)\ 



^ n =0m=0 



x sin[(n + l)$] sin [(m+1) $] sine [(n+l)e/2]sinc[(m+l)e/2]}, (29) 
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where the function f nm (x;s), x = J 7 /(y2\F\), is closely related to the sam- 
pling function 

(nlKiJ 7 , cp; -s)\m) = f nm (x; s)exp[i(n - m)cp) (30) 

for measuring the signal-mode density matrix in the photon-number basis 
11, 12]. From inspection of Eqs. ( p6|) and (|29|) we see that the symmetry 



relations 

K e ($>, il>, T, cp+n; s) = K e ($, -J 7 , V ; s), (31) 

K e ($,iJ>,r,<p+Tr,s) = K t {$-ir,il),F,y-s), (32) 

K e (-$, 1>, J 7 , V ; s) = K € ($, ip, T , V ; s), (33) 

K e {$, -V, J 7 , -<p; s) = # e ($, ip, J 7 , <p; s), (34) 

K e ($, *p, J 7 , cp; s) = K e {$, 0, J 7 , cp-iP; s) (35) 
are valid. Hence knowing the function 

K e ($,F,<p;8)=K e ($,il> = 0,r,ws), (36) 

with $, ip G (7r/2)-intervals, the function^ ($, ip, J 7 , (p; s) is known for all val- 
ues of <3>, ip, and (p. In particular, the functions K^((j), J 7 , ip; s) and K*(<fi, J 7 , ip; s), 
respectively, that are required in order to relate the the cosine and sine phase 
states distributions p c {4>, e) and p s (<f), e) to the field-strength distributions are 
given by 

K c e (<t>,r,<p;s) = K e ($ = <P,r,<p;s), (37) 
Kl^ J 7 , ip; s) = K e (<$> = 0-±7r, T, ip-\<K; s). (38) 

4 Direct sampling of the Susskind-Glogower 
phase distributions 

The function f nm (x; s) in the series expansion (p9| ) of the sampling function 
K e (&, J 7 , ip; s), Eq. fl36|), has been studied in a number of papers and different 



algorithms for numerical calculations have been discussed (see |13|, [LJ] and 
references therein). For the sake of transparency let us restrict attention to 
perfect detection (77 = 1 and hence s — 1 — 77 -1 = 0). In this case, f nm (x) = 
f nm (x, s — 0) can be written as 

fnm(x) = ^ [lp n (x)(j) m (x)) if Ul > U (39) 

[fnm{ x ) — fmn{x) if to < n), where if) n (x) and <p m (x), respectively, are the 
regular (normalizable) and irregular (unnormalizable) solutions of the energy 



eigenvalue equation of the harmonic oscillator for the nth eigenvalue [13|, [14 



The asymptotic behaviour of f nm (x) for large values of n and m can be found 
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using the semiclassical (WKB) approximation. For the argument x within 
the classical allowed region \x\ < a n = (2n + 1) 5 the function f nm (x) (m > n) 



becomes [114] 



fnm(x) ~ - (PnPm) * Pm COS (S n + COs(S m + \irj 

-PnWn[S n +\irj sin(S m + | 



where 



Pn(x) = \2n + 1 — x 2 ) 
denotes the classical momentum and 



rX 

S n (x) = / dx' 



Pn{X 



(40) 
(41) 

(42) 



is the time-independent part of the classical action. 

For proving the convergence of the series expansion of K($, J 7 , (p; s = 0), 
it is sufficient to substitute in Eq. ( |29| ) for f nm (x; s — 0) the semiclassical 
expression (flOf). For x/a n — > the functions p n (x) and S n (x) in Eq. (|40|) 
behave like n 1//2 and — (n + |)7r/2, respectively, so that (p m (x) / p n (x))z ~ 
(rn/n)'* and cos[5' n (x) + 7r/4] ~ cos(n7r/2), sin[5' n (x) + 7r/4] ~ — sin(n7r/2). 
Hence for any e > the series expansion is expected to exist, because of 
the factor (nm)" 1 that arises from the sinc-functions. Results of numerical 
calculations of i^($, J 7 , (p; s = 0) are shown in Fig. |3|. It should be noted that 
K(Q, J 7 , <f] s) separately depends on the two phases $ and <p, whereas the 
sampling function for the London phase state distribution only depends on 
the difference phase 0. 

In order to demonstrate the feasibility of direct sampling of the SG cosine 
and sine phase state distributions from the recorded difference count statis- 
tics in balanced homodyning, we have performed computer simulations of 
measurements of the phase-parametrized field-strength distributions piJ 7 , <p) 
on an equidistant grid of points {jFj,^}. We have assumed that in the ex- 
periments the signal mode to be detected is prepared in a squeezed vacuum 

|v&) = exp{-l[e(at) 2 -ea 2 ]}|0). (43) 

In particular, we have assumed that 10 4 x 30 events are recorded. From 
the results shown in Fig. [| we see that the sampled distributions are in good 
agreement with the theoretical predictions. Their deviations from the exact 
distributions obviously result from the middle observational level chosen. 
Nevertheless, the accuracy is seen to be sufficient in order to detect the 
typical features of the exact distributions. From Eqs. (|T2[) [together with 
Eqs.fliUP and (0)] and ( ]43| ) we easily see that the cosine and sine phase 
state distributions of the ordinary vacuum (£ = 0) are given by |(0| cos0)| 2 



(2/7r)sin 2 and |(0| sin0)| 2 = (27T)" 1 cos 2 . 



respectively. In this case the 
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cosine phase state distribution has a broad maximum at = 7r/2, whereas 
the sine phase state distribution has a broad maximum at = 0. When 
the value of £ is increased (£ > 0), the "vacuum noise circle" centred at 
the origin of co-ordinates in the (complex) a phase space is squeezed to an 
ellipse with the small semi-axis parallel to the real axis [[HJ. Hence, the 



cosine phase distribution is expected to become more sharply peaked at (j> 
= 7r/2. Accordingly, the sine phase state distribution is expected to show a 
double-peak structure, with the maxima close to = ±7r/2 and the minimum 
at 6 = 0. 



5 Summary and conclusions 

In the paper we have studied the problem of measuring the SG cosine and sine 
phase state distributions in balanced homodyne detection. We have shown 
that they can be directly sampled, in a coarse-graining approximation, from 
the recorded difference- count statistics with sufficiently well accuracy. The 
accuracy is determined by the interval of field strengths used for "probing" 
the phase statistics of the signal-mode under consideration. With decreas- 
ing coarse-graining parameter, this interval is increased and the accuracy is 
increased as well. 

It is worth noting that the SG cosine and sine phase state distributions 
can be regarded as special cases of ^-parametrized phase state distributions. 
The latter are based on -^-parametrized phase states that for ip = and ip = 
|7T reduce to the cosine and sine phase states, respectively. Accordingly, the 
sampling functions for the cosine and sine phase state distributions can be 
obtained by specifying the sampling function for the ^-parametrized phase 
distributions. Their integral relation to the field-strength distributions reveal 
that the sampling function exhibits a number of symmetry properties that 
can advantageously be used in calculations. 

We have calculated the sampling function using an expansion in terms of 
the matrix elements of the corresponding operator integral kernel in the Fock 
basis. These matrix elements, which are the sampling functions required for 
measuring the signal-mode density matrix in the Fock basis, can be calculated 



applying the algorithm developed in [14]. In this way, the sampling function 



can be obtained for any (non-zero) coarse-graining parameter. 

In order to give an example, we have performed computer simulations of 
measurements of the cosine and sine phase state distributions at a middle 
observational level (i.e., for a not extremely small coarse-graining parameter 
and a realistic grid of sampling points). Assuming that the signal mode is 
prepared in a squeezed vacuum state, the measured cosine and sine phase 
state distributions have been found to show the typical features of the exact 
distributions. Apart from the maxima and minima indicating the regions of 
nonclassical noise reduction and enhancement, respectively, qualitatively dif- 
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ferent shapings of the two distributions are observed which also demonstrate 
the quantum character of the state of the signal mode. It should be pointed 
out that the different shapings reflect the fact that for quantum states the 
(noncommuting) quantities cosine and sine cannot be related to a common 
phase. This effect is of course not observed in the London state phase dis- 
tribution. Clearly, the introduction of a phase that cannot be related to a 
Hermitian operator is the price that is to be paid for a common phase. 
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Figure 1: The coarse-grained cosine phase state distributions of a mode 
prepared in a coherent state \a) with (n) — 1 (a = 1) (a) and (n) = 2 (a 
= \/2) (b) are shown for e = 0.4 (dashed lines) and e = 0.8 (dotted lines). For 
comparison, the exact distributions that are observed in the limit e — > are 
also shown (solid lines). 
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(a) 
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Figure 2: The coarse-grained sine phase state distributions of a mode pre- 
pared in a coherent state \a) with (n) = 1 (a = 1) (a) and (n) —2 (a — y/2) 
(b) are shown for e = 0.4 (dashed lines) and e = 0.8 (dotted lines). For com- 
parison, the exact distributions that are observed in the limit e — > are also 
shown (solid lines). 
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(c) (d) 




Figure 3: The dependence on T and tp of the sampling function K e (<&, J 7 , <p; s) 
is shown for various values of $, the values of s and e being s = (perfect 
detection) and e = 0.4; $ = ±7r (a), $ = ±7r (b), $ = §7r (c), $ = |7r (d). 
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Figure 4: The measured Susskind-Glogower cosine (a) and sine (b) phase 
state distributions (double dotted lines) of a signal mode prepared in a 
squeezed vacuum state, Eq. (|4~3|) , with mean photon number (n) — 1 (£ = 
0.88) are compared with the calculated distributions (dashed lines), the val- 
ues of s and e being s = (perfect detection) and e = 0.4. The solid lines 
represent the exact distributions that are observed in the limit e^0. In the 
computer simulation of measurements 10 4 x 30 events (30 phase values) are 
assumed to be recorded. 



15 





